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m ■ Abstract 

I The phase space shcing method of two cutoffs for next-to-leading-order Monte-Carlo style QCD 

CN ■ corrections has been applied to many physics processes. The method is intuitive, simple to im- 

plement, and relies on a minimum of process dependent information. Although results for specific 
. applications exist in the literature, there is not a full and detailed description of the method. Herein 

such a description is provided, along with illustrative examples; details, which have not previously 
been published, are included so that the method may be applied to additional hard scattering 



p • processes. 
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I. INTRODUCTION 



Perturbative quantum chromo dynamic (QCD) calculations are essential in the effort to 
describe large momentum transfer hadronic scattering processes. At one time it was sufficient 
to work at lowest order for the hard scattering subprocesses and utilize the leading-logarithm 
approximation to treat the higher order gluon radiation and quark- ant iquark pair produc- 
tion which give rise to the scale dependence of the parton distribution and fragmentation 
functions, and to the running of the strong coupling as- As the experimental systematic 
and statistical errors decreased, the need for increased precision for the theoretical calcula- 
tions became apparent, leading to the widespread use of next-to-leading-order expressions 
for the hard scattering subprocesses with the remaining higher order terms being treated 
in the next-to-leading-logarithm approximation. Early calculations of this type were typi- 
cally performed with a combination of analytic and numerical integration techniques. The 
phase space integrations at the parton level were often performed analytically, and the con- 
volutions with the parton distribution or fragmentation functions done numerically. This 
approach is satisfactory for fully or singly inclusive cross sections, but information is lost 
about quantities over which the integrations have been performed. Thus, if cuts are to be 
placed on two or more partons (or hadrons or jets), the calculation must be started anew. 
Furthermore, for some observables it is difficult to calculate the appropriate Jacobian for 
the transformation from partonic to hadronic variables. For these reasons it was recognized 
that Monte Carlo techniques would be useful for such calculations. The Jacobians would 
be handled by the choice of histogramming variables and several observables could be his- 
togrammed simultaneously. Additionally, it would be simple to define jets and to implement 
experimental cuts on the four-vectors of the produced partons. 

In light of the above observations, a method for performing next-to-leading-logarithm 
calculations using Monte Carlo techniques was developed Two cutoff parameters serve 
to separate the regions of phase space containing the soft and coUinear singularities from 
the non-singular regions; nowadays this is referred to as the phase-space slicing technique. 

The usefulness and generality of the method may be appreciated by considering the many 
physics processes to which it has been applied. The basic core of the method was first devel- 
oped to study QCD corrections to dihadron production [Q]. It has subsequently been applied 
to direct jet photoproduction O], hadronic photon-jet [Q, direct photon W ||^, ZZ [Q, 



WW 0, WZ H, two photon [|, Z7 [TT|], and W-Higgs ^ production, nonstan- 
dard three vector boson couplings in W7 WZ [|r3| and WW [jTB| production, hadronic 
photon-heavy quark production [0, |l8l, jet photoproduction |jl9l, quantum electrodynamic 



(QED) corrections to hadronic Z production |2^, QCD corrections to slepton pair produc- 
tion ||2l|, electroweak corrections to W production p2|, single-top-quark production p3| , 
and dihadron production |^ . 



Despite this usefulness, a full and detailed description of the method does not exist in 
the hterature. Here we provide such a description. Naturally, as the method was applied to 
the above physics processes, refinements were made. We therefore take this opportunity to 
modernize and systematize the presentation relative to that given in [|I|, and show details, 
which have not previously been published. Searches for signals of new physics often rely 
on next-to-leading-order Monte Carlo-based calculations. It is anticipated that the details 
provided here will prove to be helpful for anyone wanting to apply the method to additional 
processes. 

In the course of a next-to-leading-order calculation ultra-violet singularities show up in 
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loop integrals where the momenta go to infinity. They are removed through the process 
of renormalization. (See, for example, Refs. |25, for a discussion.) Soft (infrared) 



divergences arise if the theory includes a massless field like the photon in QED or the gluon 
in QCD. They are encountered in both loop and phase space integrals and are found in the 
low energy region where the integration momenta go to zero. The soft singularities cancel 
between the virtual and bremsstrahlung processes If the massless field couples to 

another massless field, or to itself, a collinear (mass) singularity may occur in both loop and 
phase space integrals. Final state mass singularities cancel when summed over degenerate 
(experimentally indistinguishable) final states according to the theorem of Kinoshita-Lee- 
Nauenberg 1^^. For tagged hadrons there is no final state sum, and the associated mass 
singularities are factorized into fragmentation functions. Similarly, initial state singularities 
do not cancel because there is typically no sum over degenerate states; they are removed by 



factorization |28, 29 



The goal of the practitioner of next-to-leading-order calculations is to organize the soft 
and collinear singularity cancellations described above without loss of information in terms of 
observable quantities. The phase space slicing method provides a relatively simple and robust 
method to do this. Several other methods for handling the organization of the cancellations 
exist in the literature and have been used to study a wide variety of high energy processes, 
including some of those listed above. The phase space slicing method of one cutoff first 
developed in 13^, |31|] divides the phase space according to Sij = {pi + Pj)"^ > ysi2 where 
Pi and Pj label the momenta of partons i and j, and y is a small dimensionless parameter. 
Another variant for jets |3^, Q and hadrons and heavy quarks |0 partitions phase space 
according to Sij > Smin where Smin is a small dimension-full parameter. It is also possible 
to engineer the singularity cancellation using plus distributions, commonly referred to as 
the subtraction method, which has been applied to jets ^ ^ ^ 

and heavy quark final states H3, HI, H6[. The subtraction method taken together with 



factorization formulae that interpolate between the soft and collinear approximations to the 
matrix elements is known as the dipole method The dipole method was originally 

developed for jet and light hadron cross sections where it has seen extensive application 
p8| , 49, 15^, 51, 52, 53, 54, 55, 56|. It has recently been extended to handle massive fermions 
and partons ||5^, |5^, |59| and is finding many applications |]23| , |60| , |6T1 , |62| in that domain as 
well. A brief comparison of the slicing and subtraction methods is given in Appendix A. 
Properly implemented, all methods should give identical physics predictions. 

This paper proceeds as follows. In Section II we give details of the phase space slicing 
method with two cutoffs. We examine the soft and collinear regions of phase space and 
see how to arrive at a finite cross section. In Section III the process of electron-positron 
annihilation into quarks is studied for massive, massless, and tagged final states. After that, 
the examples of lepton pair production and single particle production in hadronic collisions 
are given. We conclude in Section IV. As mentioned above. Appendix A contains a brief 
comparison with the subtraction method. Appendix B contains angular integrals useful in 
the soft analysis. A discussion of terms that vanish like the ratio of the two cutoffs is given 
in Appendix C. Appendix D explains how to improve numerical convergence. 



II. THE METHOD 

This section contains the main derivations for jet, fragmentation, and heavy quark final 
states, as well as a discussion of initial state mass factorization. Before getting into too 
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many details, it will be helpful to outline the procedure first. The typical calculation involves 
lowest order two-to-two subprocesses which have two-body final states and higher order two- 
to-three subprocesses which lead to both two- and three-body final states. In addition, the 
one-loop virtual corrections also contribute to the two-body final states. 

We begin by decomposing the three-body phase space used to calculate the two-to-three 
contribution to the partonic cross section into two regions which we call soft, S, and hard, 
H, by writing ^ 

where $ = A^/^(s,m^, ml) is the usual flux factor which depends on the partonic center-of- 
momentum energy squared s and the incident particle masses mi and m2, I]|^3p is the 
two-to-three body squared matrix element averaged (summed) over initial (final) degrees of 
freedom, and dT^ is the three-body phase space. The partitioning of phase space into S and 
H depends on a parameter 6s in a manner to be described below. Within S the double pole 
(eikonal) approximation to the matrix elements is made and then analytically integrated 
over the unobserved degrees of freedom in n space-time dimensions. The result, depending 
on the masses of the partons involved, may contain double and/or single poles in n — 4, and 
accompanying double and/or single logarithms in the soft cutoff 6s. We always work in the 
approximation where the cutoffs are small, so terms of order 6s may be neglected. Just how 
small is needed will be studied below. 

Next, if there are collinear singularities present, the hard region is further decomposed 
into collinear C, and non-collinear C, regions as follows: 

2<l>iH^' ' 2$yHC^' ' 2<l>iHC^' ' ^ ^ 

This partitioning depends on a second cutoff 6c- Within HC the leading collinear pole ap- 
proximation to the squared matrix element is made. As explained below, exact collinear 
kinematics may be used to define the integration domain of HC when 6c ^ 6s. The inte- 
grations over the unobserved degrees of freedom are performed analytically in n space-time 
dimensions giving a factorized result where single poles in n — 4, and single logarithms in 
both cutoffs 6c and 6s, multiply splitting functions and lower-order squared matrix elements. 
This is in the approximation where terms of order 6c and 6s are neglected. 



The cancellation of poles in n — 4 is based on the Kinoshita-Lee-Nauenberg theorem p7 
or mass factorization [^, |29|, depending on the situation. For experimentally degenerate 
final states, the soft and final state hard collinear singularities cancel upon addition of the 
interference of the leading order diagrams with the renormalized one-loop virtual diagrams. 
The remaining initial state collinear singularities are factorized and absorbed into the parton 
distribution functions. The result is finite in n = 4 dimensions, but depends logarithmically 
on the cutoffs. For tagged hadrons there is no final state sum, and the associated mass 
singularities are factorized into fragmentation functions. 

The integration over the hard non-collinear HC portion of the phase space is performed 
using standard Monte Carlo techniques. The result is finite by construction and the ex- 
pressions may be evaluated in four dimensions. At this stage, the calculation yields a set 



^ Implicit in the cross section is a measurement function which serves to implement the jet algorithm and/or 
define the experimentally visible portion of phase space (the cuts) . For the cancellation of singularities to 
take place, the measurement function is required to be infrared-safe pJ. 
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of two-body weights which have exphcit logarithmic dependence on the two cutoffs and the 
three-body weights for which a logarithmic dependence on the cutoffs develops as the Monte 
Carlo integration is performed. When all of the contributions are combined at the histogram- 
ming stage, the cutoff dependence cancels for suitably defined infrared-safe observables. In 
the following subsections we look at each of these steps in detail. 



A. Soft 

In this subsection we describe the procedure for extracting soft gluon singularities. When 
one of the gluons is soft, the phase space is greatly simplified and the eikonal (double pole) 
approximation of the matrix element is valid |6^, |6^, 0, The cross section is 

simple enough to be analytically integrated over the unobserved degrees of freedom in n 
space-time dimensions. The required integrals are well known P5| , |69| , 0. The result. 



depending on the masses of the partons involved, contains double and/or single poles in 
n — 4, and accompanying double and/or single logarithms in the soft cutoff 6s- We work in 
the approximation where terms of order 6s are neglected. 

Generically, we begin by writing the two-to-three body contribution to the partonic cross 
section 

^ = ^/El^3l'rfr3, (2.3) 

as the sum of soft, S, and hard, H, terms 

a = as + (Th , (2.4) 

where 

'^s = ^lTm'dT,, (2.5) 

and 

^n = ^I^T.\M3\'dT,. (2.6) 

In this section we examine as in detail. Further evaluation of cth is deferred to the following 
subsections. 

Let the particles in the scattering be labeled by their four-momenta pi+p2 = Ps +P4 +P5 
and define the Mandelstam invariants Sij = {pi + Pj)'^ and tij = {pi — pjY. Consider the 
case when parton 5 is a soft gluon. The soft region S is defined in terms of the gluon 
energy in the pi + p2 rest frame by < -E5 < 6s^/sl2/'^■ The hard region H is the 
complement: E<^ > 6sy/si2/2. The gluon energy can be calculated from the other invariants 
in the problem as follows. Start with P1+P2—P5 = P3 +Pa which, after squaring both sides, 

yields (pi+P2)^-2p5-(Pi+P2) = (^3+^4)^- In the pi+p2 rest frame pi-Fp2 = v^(l'O'O'O)' 
so s 12 — 2E5^/si2 = S34. Solving for the gluon energy gives 

E,= '-^§^. (2.7) 



This expression for E^ and the definition of the soft region are independent of the masses 
of the other particles in the reaction. Now that we have defined the boundaries of the soft 
portion of phase space, we examine the approximations that can be made. 
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The three-body phase space in n dimensions is given by 
d"-^P3 dJ'-^pi dJ'-^pr, 



2p§(27r)"-i 2p^(27r)"-i 2pg(27r) 



^(27r)V(pi+P2-P3-P4-P5) 



(2.8) 



The divergence in the integral of the matrix element over phase space will be at worst 
logarithmic. Therefore, up to corrections of 0{5s)., we can set = in the delta function 
and regroup terms, yielding 



c^rslsoft - 

This is simply 
where 



2p^(27r)"-i 2pl{2nY 



-(27r)V(pi+p2-P3-P4) 



n-1 



dV 



3 1 soft 



dV. 



d'^-^p^ 



2pI{2tiY-^ ' 



dV, 



d^ d"- Pi 



f3^(27r)"5"(pi+p2-P3-P4), 



(2.9) 



(2.10) 



(2.11) 



2p§(27r)'*-i 2p^(27r)' 

is the two-body phase space of partons 3 and 4. Likewise, up to corrections of 0{5s), we 
can parameterize the gluon's n— momentum in the pi + p2 rest frame as 



P5 = -£'5(1, . . . ,sin^i sin ^2, sin ^1 cos ^2, cos ^1) , 



(2.12) 



where the dots indicate the n — 4 unspecified momentum components which may be trivially 
integrated over using 



d^-'p^ = dmp^\^-Mn^_2 



dE^E^-^ sin"-=^^i dOi sin"-^^2 rf^2^^n-4 ■ 



n-4/ 



The angular volume element 



^n— 4 — 



may be rewritten using 



r[(n - 3)/2] = V7r2 



r[(n - 3)/2] ' 

rfl-2e) 



2e' 



(2.13) 
(2.14) 
(2.15) 



r(l - 6) ' 

where we have set n = 4 — 2e. The final result for the phase space volume approximated in 
the soft region is 



dr 



3 1 soft 



dVo 



Ativ r 1 



5127 r(l - 2e) 2(27r) 



dS, 



with 



dS 



' ^ ' dE5El-^'sm^-'^'9id9ism-^'92d92. 



= -{-) I 



(2.16) 



(2.17) 



Once we have the corresponding soft approximation to the matrix element this integral can 

be performed, yielding the advertised singularity structure. But before proceeding, we pause 
to note that choosing 5 to be the soft gluon is not special. A similar analysis holds when 
parton 3 or 4 is a gluon. At this order in perturbation theory, only one gluon may be soft 
at a time. 
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Soft photon emission in QED is characterized by the factorization of an eikonal current 
from the scattering amphtude. The structure of multiple soft photon emission has been 
studied by Grammer and Yennie |Q. In QCD the process is different because gluons 
carry color charge and can therefore radiate during the scattering. Fortunately, when QCD 
amplitudes are decomposed into color sub- amplitudes they enjoy the same factorization 



properties as QED amplitudes |B1 

Let parton 5 be the soft gluon and take it to have color index a (= 
Lorentz index yu. The matrix element factorizes as 



3 I soft 



_ 1) and 



(2.18) 



The mass dimensions of the strong coupling have been isolated into the parameter which 
is identified with the renormalization scale, leaving the dimensionless coupling g. The soft 
gluon's polarization vector, denoted by ^^(ps), is Lorentz contracted with the non-abelian 
eikonal current given by 

4 „M 



(2.19) 



which itself is color contracted with the color sub- amplitude M2. The sum corresponds 
to the soft gluon being emitted from each external line in turn. The SU{N) color charge 
associated with the emitting parton / is denoted by Tj. Squaring and summing Eq. ( p.l8| ) 
over the soft gluon polarizations gives 



where 



soft 



2 2e 
-9 f^r 



Pf ■ Pf 



fjt^ PfP5 Pf -PS 



(2.20) 



(T°.M2)(T'^,M2 



Cl. 



b^/ ...C4 



rria rpa 



M, 



Cl. 



i-f ...d^'l ...C4 ; 



(2.21) 



is the square of the color connected Born amplitude. If the emitting parton is a final 
state quark or initial state anti-quark, the color charge is in the fundamental representation 



, A^). For a final state anti-quark or initial state quark (T"'] 



If the emitting parton is a gluon, the color charge is in the adjoint representation {T" 
Substituting (gl|) and into (g3|) gives 



Ffl- 





27rF(l-2e) 1^ 



-Pf ■ Pf 



PfP5 Pf -Pb 



-dS . 



where 



1 



— yM° ,dF2 



(2.22) 



(2.23) 



The integration over the eikonal factors depends on the masses of the particles in the reaction. 
We leave the integrals to be performed on a case-by-case basis, although all possible mass 
combinations may be worked out. Specific applications of Eq. ( p.22|) are given below. 

From a practical point of view, the presentation just made requires only the calculation 
of the eikonal factors and the color connected Born amplitudes in n = 4 — 2e dimensions. 
There is no need to evaluate the full two-to-three body matrix element in n dimensions. 
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Depending on the complexity of the process, this may be a simphfication. However, if the 
full two-to-three body matrix element in n = 4 — 2e dimensions is at hand, setting = 
everywhere in the numerator and retaining only the leading singular terms as — will 
reproduce Eq. ( |2.20|) directly, this being known as the double pole approximation 



B. Collinear 

We now return to the further evaluation of the hard portion of the cross section which 



was separated out in Sec. |II A] . The phase space is greatly simplified in the limit where two 
of the partons are collinear. In the same limit, the leading pole approximation of the matrix 
element is valid. The cross section is simple enough to be analytically integrated over the 
unobserved degrees of freedom in n space-time dimensions. The result contains single poles 
in n — 4, and accompanying logarithms of the soft 5s and collinear 5c cutoffs. Terms of order 
5s and 5c are neglected. 

To this end we further decompose (Th given in Eq. ( p.6| ) into a sum of hard-coUinear HC 
and hard-non-coUinear HC terms 



^HC = ;i/ ElMal'cirs, (2.25) 



<yn = 0-Hc + cthc ' (2.24) 

with 

2$ 7hc^ 
and 

^HC = 4^El^3rrfr3. (2.26) 

The HC regions of phase space are those where any invariant (sjj or tij) becomes smaller 
in magnitude than 5cS\2^ the collinear condition, while at the same time all gluons remain 
hard. The complementary HC pieces are finite and may be evaluated numerically in four 
dimensions using standard Monte-Carlo techniques ||71 . 



The piece containing the collinear singularities, cthc; is treated according to whether 
the singularities are initial or final state in origin. For the former, factorization provides 
the formalism for removing the singularities. In the latter case, we distinguish between 
experimentally degenerate and tagged final states, and rely on either the Kinoshita-Lee- 
Nauenberg theorem or factorization to dispose of the singularities. We discuss the final 
state cases first, then return to the initial state. 



1. Indistinguishable Final States 



Consider the case when there is a sum over experimentally degenerate final states, such 
as a jet or total cross section. Let partons 4 and 5 be massless and collinear to each other, 
< S45 < 5cSi2- If we define P45 = + p^, then for fixed ps we have d^'^p^^ = d^'^p^. The 



three body phase space Eq. ( |2.8| ) may be written as 



dV 



3 coll 



2p[](27r)"-i 2p%{2TiY 



-(27r)V(pi+p2-P3-P45) 



d 



n-l 



2pg(27r)"-i pI 



(2.27) 
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This is simply 



(iralcou - Q, . , , (2.28) 

where dT2 is the two-body phase space of the particles 3 and 45. In the collinear limit 
{pt —>■ with z fixed) we can write 



P45 
P5 



(P,0,0,P) 



izP + ^,Pt,zP) 
iil-z)P + 



Pt 



2(1 -z)P 



(2.29) 



Then P4 + P5 = P45 + ^iPt) and 



S45 = 2p4 ■ P5 



Pt 



z(l-z) 



(2.30) 



Now using d"" p = dpp^ diln-s and z{l — z)ds45 = dpf we find 



^" P45 
2pg(27r)"-i p2 



(4vr) 



levr^rri - e 



-(i^; (is 45 [545^(1 — -z)]' 



(2.31) 



The corresponding approximation to the matrix element is obtained by imposing collinear 
kinematics on the portion of the two-to-three matrix element proportional to the leading 
collinear singularity. This is known as the leading pole or collinear approximation. As a 
consequence of the factorization theorems the squared matrix element factors into 



the product of a splitting kernel and a leading order squared matrix element [0, |72 
As above, let partons 4 and 5 be collinear. The matrix element factorizes as 



5^|M3(1 + 2 ^ 3 + 4 + 5)p ^ ^|M2(1 + 2^3 + 4')\^P,,iz, e)g'fi. 



2,.2e_ 

r 



S45 



(2.32) 



where the Pij{z,e) are the unregulated {z < 1) splitting functions calculated in n = 4 — 2e 



dimensions related to the usual Altarelli-Parisi splitting kernels [|72[. We label as 4' the 
parton which splits into the 45 collinear pair. Generally, Eq. (|2.32| ) contains an additional 
term that vanishes after integration over the azimuthal angles in n dimensions ||3^, ^ . Such 
a term does not contribute to our result. 

Substituting Eqs. (g^), (|3l|), and (|23^ ) into Eq. (121251) gives 



j-.l+2^3+4+5 

acr jjc 



a, rri 



X 



27rr(l 

S12 



2e) 



S12 



■S45 



J dz[z{l - z)]-'P^i,{z, e) , (2.33) 



where we have used 



1 



r(i-e) 



(2.34) 



r(l-e) r(l-2e) 

The collinear condition (0 < S45 < ScSi2) sets the S45 integration limits. The hard condition 
sets the z integration limits which depend additionally on the splitting function involved. 
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and also on the mass of the parton recoihng against the 45 pair (3 in this case). This later 
dependence enters through the threshold condition, as discussed below. 

First take parton 3 to be massless. For q ^ qg splitting the hard region is defined by 
^s^/svi/'^ < E's < ^^12/2, assuming 4 labels q and 5 labels g. From Eq. ( p.29| ) we have 
■S34 = (P3 + Pif = 2^3 ■ Pa - (2^3 ■ P45)z and S12 = (ps + ^45)^ = 545 + 2^3 ■ ^45 - 2p3 ■ P45 
which together yield S34 ~ zsi2. Using Eq. ( |2.7| ) the hard condition becomes 

0<z<l-5s. (2.35) 

For the g ^ gg splitting it is required that both gluons be hard, i.e., E4 and > Ss^/svi/'^- 
z then satisfies the relation 5s < z < 1 — 5s. For the g ^ qq splitting there are no soft 
singularities, so < 2; < 1 may be taken. In all of these cases the z integration limits 
are independent of S45 by virtue of the approximation S45 = implicit in S12 — 2p3 • P45 
(this point is discussed further at the end of Sec. |111 B| and in Appendix C). The outermost 
integration over S45 may therefore be performed giving the result 



J 1+2^3+4+5 



as 



Til 



27rr 1 -2e 



'47r/4 



- 1 5: 



dzz 



For z < 1 the n— dimensional unregulated splitting functions may be written as Pij 
P.,,{z) + ePl^{z). Explicitly, 

■ Up 



Pw{z,e) . 

(2.36) 

z,e) = 



P,9i^ 



1 - z 
-Cpil - z) 



C, 



1 + (1 



-Cpz 



2N 



1 - z 



+ 



+ z(l-z) 




1 
2 



K + (1 - ^; 

-z{l - z) , 

where = 3 and Cp = {N'^ — l)/2N = 4/3 for QCD. Expanding the integrand in Eq. 
to 0{e) and integrating over z yields the final state hard-collinear terms 



^-.l+2-.3+4+5 



a. 



r(l-e) [AtthI 



271 T(l -2e) 



S12 



l4'-»45 



/14' 

^ — + Af!-^' 



where 



At 

At 
A{- 



+99 

*qq 
*99 
*99 



Cj. (3/2 + 2 In 5, 



c 



7/2-7rV3-ln2 5, - In 4 (3/2 + 2 In 5, 

n//3 (In 5c -5/3) 
(11/6 + 2 In 5,) 



A^ 



67/18 - 7rV3 - In^ 5, - In 4 (11/6 + 2 In 5s) 



2.37) 
2.38) 

2.39) 
2.40) 
2.41) 
2.42) 

2.43) 
2.44) 

Pi 



2.45) 

2.46) 
2.47) 
2.48) 
2.49) 
2.50) 
2.51) 
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where n/ denotes the number of active flavors. 

When the mass of the parton recoihng against the 45 pair, 1713, is retained, the factoriza- 
tion of the phase space and matrix element is unaffected. Likewise, the colhnear condition 
< S45 < ScSi2 remains unchanged. It is only the boundries of the hard region that are 
modifled. The full kinematic range of the invariant S34 is < S34 < S12. The thresh- 
old for producing a particle of mass sets the lower limit. In terms of S34 the hard 
region is < S34 < {1 — 6s)si2. This implies the hard condition Eq. ( p. 351 ) becomes 
< z < 1 — Ss/{1 — m|/si2). A similar analysis follows for the g ^ gg splitting case 
wherein 6s/ (1 — m|/ S12) < z < 1 — 6s/ (1 — ml/ S12). We therefore immediately see that Eqs. 
( p.46| )- (|2.51| ) are valid with the replacement 6s —>■ 6s/{l — rnf/su)- 



2. Tagged Final States 

Next, consider a process where a particular type of hadron is identifled in the flnal 
state. This necessitates the introduction of a fragmentation function Dh/c{z) which gives 
the probability density for flnding a hadron h which carries a fraction z of the momentum 
of the parent parton c. Consider the case where parton 4 fragments into a hadron h, for 
which the lowest order cross section is 



Dh/A{z)dz. 

The hard collinear cross section expression in Eq. ( p.36|) becomes 



(2.52) 



as r(i 



'47r/i^ 



6: 



X 



dyy~ 



27rr(l -2e) ^ S12 / 
yr^Dh/4{x)dxP44'{y,e)6{xy - z)dz . 



(2.53) 



The delta function insures that the hadron h carries a momentum fraction z of the parent 
parton's momentum (parton 4' in this example). Here there is a splitting 4' — > 45 followed 
by parton 4 fragmenting to hadron h. When all possible 2 — > 3 subprocesses are considered, 
there will be several contributions of this same form, corresponding to a sum over the parton 
4. For example, if 4' is a gluon, there can he g ^ gg followed hj g h ot g qq followed 
hj q ^ h or q ^ h. Similarly, if 4' is a quark q, there can he q —* qg followed by g — /i or 
hy g ^ h. Furthermore, the limits of integration on y depend on the splitting function as 
in the case discussed in the previous subsection. 

The collinear singularity, evidenced by the pole in e in Eq. (|2.53|) , must be factorized and 
absorbed into the bare fragmentation function. To do this, we introduce a scale dependent 
parton fragmentation function 



Dh/c{z,Mj] 



D 



h/c 



as r(l - e) / 47r/i2 
2^r(l -2e) 



dy 

y 



Dn/c'{z/y)PX{y). (2.54) 



In this expression there is an implied sum over the index d corresponding to the sum over 
the different fragmentation possibilities referred to above. The flnal state factorization scale 
has been denoted by Mj. Notice, too, that the integration over y extends from 2; to 1. This 
form for the scale dependent fragmentation function corresponds to the MS convention. The 



11 



regulated (x < 1) splitting functions [^] are given by 
Cf 



1 + (1 -x)^ 



X 



2N 



X 



'l-x) 



H h xil — x) 



X 



+ C-lN-\nf]5{l-x) 



1 r 
2 



(2.55) 
(2.56) 
(2.57) 
(2.58) 



Next, we rewrite the bare fragmentation function in Eq. (|2.52|) in terms of the scale dependent 
expression given above, yielding to 0{as) 



X {Dh/iiz,Mf)+ { - 



as r(l-e) /47r^^ 
2^r(l-2e) V Mj ^ 



(2.59) 



The second term is sometimes refered to as the mass factorization counterterm. When dao 
and dauc are added together, there is a cancellation between the two singular expressions. 
Note, however, that this cancellation is not complete since the limits of the y integration in 
the two expressions differ. 

After the cancellation, the resulting 0{as) expression for the fragmentation contribution 

is 



OLs r(l — e) ( Airfxl 



X \Dh/4'{z,Mf) + 



2ttT{1 - 2e) V si2 / 
Af'^i' ^4 + 5) 



+ Ag^(4' ^4 + 5) 



Dh/A{z,Mf)jdz. 

(2.60) 



The soft coUinear factors A^'^ result from the mismatch in the y integrations in the fragmen- 
tation and subtraction pieces, mentioned above. They are given by 



ASC 

^0 



ATiq-^qg) = C^(21n55 + 3/2) 
ATig^gg) = 2N\n6s + {UN - 2nf)/6 . 



The modified fragmentation function Dii/c{z, Mf) is given by 
where 



-jD,/,{z/y,Mf)P^,:^. 



P,:' = Pc'ciy) In 



y{l - y)6cSi2 
M) 



P'c'ciy) ■ 



(2.61) 

(2.62) 
(2.63) 



(2.64) 



(2.65) 
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P{y) and P'{y) are the n = 4 and (9(e) pieces, respectively, of the unregulated splitting 
kernels given in Eqs. (|2.37|) - (|2.44|) . The D functions contain an explicit logarithm of 6c as 
well as logarithmic dependences on 6s which are built up by the integration on y when c' = c. 
In Appendix D it is shown how to make the Sg dependence explicit, thereby improving the 
convergence of the Monte Carlo integration. 

Comparing with the previous subsection, we see that when going to the fragmentation 
case, the hard collinear terms, Eqs. (|2.46|) - (|2.51|) , for the fragmenting parton are replaced 



by a combination of the D function and soft collinear factors A^'^. Nevertheless, a careful 
comparison of the two cases shows that the poles in e cancel and the final results for physics 
observables are independent of the cutoffs. This will be illustrated by several examples to 
follow. 



3. Initial State 

The treatment of the initial state collinear singularities is much the same as that for the 
previous case of final state fragmentation. The collinear singularities are absorbed into the 
bare parton distribution functions leaving a finite remainder which is written in terms of 
modified parton distribution functions. In addition, there are accompanying soft collinear 
factors as in the fragmentation case. However, some of the details are different, so a brief 
summary of the derivation is given here. In order not to unnecessarily complicate the 
discussion, only the details for one of the incoming partons will be shown. 

Consider a process which involves a parton on leg 2 coming from an incoming hadron B, 
so that in lowest order 

^^0^+^^=^+^ = G2/B{x)dxdal+^^''+'' , (2.66) 

where G2/B{x)dx denotes the probability of getting parton 2 from hadron B with a momen- 
tum fraction between x and x + dx. The hat symbol" is used here to label a purely partonic 
subprocess. We are interested in the next-to-leading-order corrections coming from the vari- 
ous possible parton splittings which can occur on leg 2. The hard collinear contribution Eq. 
( p.25| ) is calculated by applying the collinear approximation to the appropriate three-body 
matrix elements as follows: 

Y}M^{1 + 2 ^ 3 + 4 + 5)p ~ El^2(l + 2' 3 + 4)|^P2'2(^, e)^?^^^ , (2.67) 

where z denotes the fraction of parton 2's momentum carried by parton 2' with parton 5 
taking a fraction (1 — z). Using the approximation P2 — P5 — zp2, the three-body phase 
space may be written as 



dT 



3 coll 



2p0(2vr)"-i2p0(2vr)«-i('"^ ^ (pi + - Pa - P4) 



(2.68) 



in— 1 



The square bracketed portion is just two-body phase space evaluated at a squared parton- 
parton energy of zsi2. The dependent part may be rewritten as 

^ ^ -dzdt2,[-{l-z)t2,r. (2.69) 



2p^(27r)"-i 167r2r(l - e) 
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The allowed range for ^25 is given by the collinear condition < —^25 < ^cS 12 ■ The ^25 
integration yields 

1 



-dt 



--(5cSl2) 



(2.70) 



Using these results, the three-body cross section in the hard collinear region may be written 

as 



G2/B{y)dydao 



l+2'^3+4. 



ZS12, tl3, tu) 



27rr(l -2e) \ s 



12 



X 



- ) 5^ ''P2'2{z-, e)dz{l — z) "dd/z — x)dx . 



(2.71) 



Note that a factor of 1/z has been absorbed into the flux factor for the two-body subprocess. 
The delta function insures that the fraction of hadron B' s momentum carried by parton 2' 
into the two-body subprocess is x in order to be able to combine this result with the lowest 
order contribution. The delta function may be used to perform the y integration, but one 
point must first be made. S12 is related to the square of the overall hadronic squared center- 
of-mass energy S by S12 = yS. On the other hand, in the lowest order subprocess the relation 
is S12 = xS. It is convenient to rewrite the above expression using this latter definition for 
S12. Therefore, after the y integration, we obtain 



,1+5^3+4+5 



G2/B{x/z)dal'^^'^^^^{si2, tl3, tu) 



as r(l - e) ( A-Kfi. 



X 



SrP2'2{z,e 



dz 

z 



27rr(l - 2e) 
dx . 



S12 



(2.72) 



Comparing with the corresponding result for final state fragmentation, we see that a factor 
of [z{l — z)]'"" has been changed to [(1 — ^)/^]~^ 

In order to factorize the collinear singularity into the parton distribution function, we 
introduce a scale dependent parton distribution function using the MS convention: 



Gt/B{x,Hf) = Gb/B{x) + — 



a, r(l-e) Unfil 
2^r(l -2e) \ 



1 dz 

—Pw{z)Gy,B{x/z). (2.73) 

z Z 



Next, using this definition, we replace G2/b{x) in the lowest order expression (|2.66| ) and com- 
bine the result with the hard collinear contribution ( p.72| ). The resulting 0{as) expression 
for the initial state collinear contribution is 



^ l+B^3+4+5 



coll 



'4:iTfi: 



X <G2'/B{z,fif) + 



2nT{l - 2e) ^ S12 / 
Af{2 2' + 5) 



+ A^^(2 2' + 5) 



G2/B(^,/i/)Hfc74) 



Note that in this expression the soft collinear factors (given in Eqs. ( |2.61| )- (p.63| )) depend 
on the initial state factorization scale fif. The G functions are given by 



Gc/B{x,flf) = 



l-SsS, 



dy — 
—G^>/B{x/y,i2f)Pcc'{y) 



(2.75) 
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FIG. 1: Leading order contribution to electron-positron annihilation via photon exchange. 



The n = 4 and 0{e) pieces of the unregulated splitting kernels, P{y) and P'{y), are given in 
Eqs. ( p.37|) -( ^.44D . An example of a hadron-hadron process will be given in the next section. 

As in the final state hadron case, the G functions contain an explicit logarithm of Sc 
as well as logarithmic dependences on 6s which are built up by the integration on y. In 
Appendix D it is shown how to make the 6s dependence explicit, thereby improving the 
convergence of the Monte Carlo integration. 

III. EXAMPLES 

In this section we provide five illustrative examples applying the method developed in the 
previous section. The results are shown to be in complete agreement with those available 
in the literature. We begin by calculating the QCD corrections to electron-positron annihi- 
lation into a massive quark pair. The quark mass serves to regulate any would be collinear 
singularities. There are only final state soft singularities and, hence, only the soft cutoff is 
required. Next, the QCD corrections to electron-positron annihilation into a massless quark 
pair are considered. In this case final state soft and collinear singularities are encountered, 
necessitating the use of both soft and collinear cutoffs. The example of inclusive photon 
production in hadronic final states of electron-positron annihilation is then presented, illus- 
trating the use of fragmentation functions. Finally, we close the examples section by showing 
how to calculate the QCD corrections to lepton pair and single particle inclusive production 
in hadron-hadron collisions. Both examples contain initial state soft and collinear singular- 
ities, necessitating the use of scale dependent parton distribution functions. Furthermore, 
the single particle inclusive cross section calculation also requires the use of scale dependent 
fragmentation functions. 

A. Electron-positron annihilation to massive quark pair 

Electron-positron annihilation into a massive quark pair has a particularly simple singu- 
larity structure, that of soft singularities in the final state only. It will therefore be used as 
a first example of the method described in the section above. 

Working in the single photon exchange approximation, the leading order Feynman dia- 
gram is shown in Fig. |I]. Neglecting the electron mass and denoting the quark mass by m. 



with 




(2.76) 
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(a) (b) 
FIG. 2: Real emission contribution to electron-positron annihilation. 

the leading order cross section 

rfa° = -^J|M2|2rfr2, (3.1) 

calculated in n = 4 — 2e dimensions is expressed in terms of the (summed and averaged) 
matrix element squared 

J2W = 2Ne'Qi n' +^-6 , (3.2) 

V ^12 Sl2 / 

and the two-body phase space 

dT, = ±- ( ^ ) P'-'-j^T^, / sin^-^^e d9 . (3.3) 
167r Vsi2/ r(l - e) Jo 



The center-of-mass scattering angle is denoted by 6 and f3 = \Jl — Am? / Si2. Qg is the quark 
charge in units of e and = 3 is the number of colors. When masses are present, it is often 
convenient to define primed Mandelstam invariants which are the ones defined previously 
minus some combination of squared masses. In this case, s[j = Sij — rn? and t[j = tij — rn?. 
Performing the integration over 9 we obtain the well known e = result, 

.o = ip!oJiV(l + H!^),3. (3.4) 

3Si2 Si2 

with a = e^/AiT. 

Because the quark mass regulates any would-be final state collinear singularity, the ap- 
propriate decomposition of the two-to-three contribution to the cross section is simply given 
by Eq. ( |2.4| ). For the 0{as) QCD corrections we therefore need the soft cross section (p. 51), 
the hard cross section ( ^.6[ ), and the virtual corrections. 

The real emission diagrams that give IM3P are shown in Fig. ^. If we define the (summed 
and averaged) squared matrix element as 

Y^\Ms\' = ]e'Qy^NCp^, (3.5) 

then 

^ = tr[^27''^i7l 

X {tr [{^3 + m)7^(^4 + ^5 - - m)7^(-^4 - ^5 + m)'y^] / 545^ 



'45 



+2 tr [(^3 + m)7''(^3 + ^5 + m)7''(^4 - m)7^(|*4 + ^5 - m)7^] / S35S4 
+ tr [(^3 + m)7"(|(3 +^5 + m)Y{-^i + m)-f,{^3 + ^5 + m)-/^] , (3.6) 
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where the strong couphng is denoted by g and Cp = {N^ — l)/2N = 4/3. To obtain the 
hard contribution, cth, the traces may be evaluated in four space-time dimensions. 

There is a soft singularity when the energy of the gluon in Fig. ^ goes to zero. The 
corresponding soft contribution to the cross section, as, is given by Eq. ( p.22| ). In this 
case, the sum in ( p^.22| ) is taken over the final state quark legs (labeled 3 and 4) and the 
color connected Born cross sections are related to the leading order cross section by dcr^^ = 
Cpda'^. We find 



daa = da^ 



a, rri 



27rr 1 - 2e 



47r/i2 

. ^12 



m 



m 



2m^ 



{P3 ■ Pf>Y {Pa ■ PbY P3 -PbPi- P5 



dS 



(3.7) 

The poles need to be integrated over the soft phase space according to Eq. ( |2.17| ) to extract 
the singularities in dimensional regularization. Define 



lis'. 



35) 



lis'. 



45) 



^35 



■^45 



dS 



dS 



■535 ■^45 



dS. 



(3.8) 
(3.9) 
(3.10) 



In terms of the pip2 center-of-momentum scattering angle 6, ps and p^ may be written as 



P3 
Pi 



2 

S12 



(1,0, ...,O,/5sin^,/5cos0) 
1, 0, . . . , 0, —[3 sin 6', —[3 cos ( 



(3.11) 



I 

^35 



2P3 ■ P5 
2p4 ■ P5 



Using these together with Eq. ( p. 12 ) we find 

/ii2£'5(l — (3 sin 6 sin 61 cos 02—13 cos 6 cos 9i] 
S45 = 2p4 ■ p5 = ^/sl2E^{l + /5 sin 6* sin 9i cos O2 + (3 cos 9 cos ^^i) 

The gluon energy integrals in Eqs. (PD-(Fl0D may be performed trivially 

1 

S12 



(-) 

VS12/ 



dE.E}-^' 



S12EI 



1 

2^ 



-2e 



(3.12) 



(3.13) 



The remaining angular integrals are well know and are tabulated in the Appendix B. The 
complete results are 



HS'S5S'J 



2e 1-/3 



1 



2m? 



h In 5, In — 

2e " 2(3 I -(3, 



Li2- 



2(3 



4 



+ In 5s In 



l + (3' 



si2(3\ 2e I -(3 I + (3 4 1-/3 I -(3^ 

We may therefore write the final expression for the soft cross section as 



cfcTc = da^ 



a. Til - e) 



'47r/i^ 



27rr(l-2e) ^ S12 



(3.14) 
(3.15) 

(3.16) 
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(c) 



(d) 



FIG. 3: Loop and counterterm corrections to electron-positron annihilation via photon exchange. 



where 



2(3 



ACf 

1 + 13' 
2(3 



In (5o H In — -^-^ 

2(3 1-/3 



Li2 



2(3 



1, 
+ - In 



l + (3 



l + (3 4 1-/5 



In 5o In 



l + (3' 



(3.17) 



The virtual contribution is obtained from the one-loop diagrams shown in Fig. ^. In the 
on-shell renormalization scheme diagrams (c) and {d) cancel exactly. The vertex correction 
needed in Fig. |^(a) is shown separately in Fig. |^. After performing the loop integrals the 
result for the vertex valid for > Am? is 



= {-leQq) 6ij (g CpCA u{p2) 



A 



u{pi 



with 



A 
B 



l + (3 



(3.18) 



(3.19) 



+ 



1+^ 

/5 



--In^ 



l + (3 



+ 21n 



1-/3^^ 2(3 
l + (3 ^l + (3 



+ 2Li2 



.1-/3 277 



1 + /3 



+ 



and 



a 



2-e 



(27r; 



4-2e 



r(i + e) 



m? 



(3.20) 



(3.21) 



18 



FIG. 4: The vertex correction. 



The vertex counterterm, Zi, imphcit in Fig. is fixed in the on-shell renormaUzation 
scheme by the condition that the renormaUzed vertex through 0{g'^) evaluated at zero 
momentum transfer equals the leading contribution {—ieA^{q = 0) = —iej'^). This results 
in 

(3.22) 



Z^ = l-g'CFa{- + 4 



The interference of the leading order diagram with the one loop renormalized diagrams 
yields 



with 



as r(l - e) / 47r/i; 
2^r(l-2e) 



^ + Al]+da 



(3.23) 



A"" 
^0 



2(3 



+ 



2fl-l±^lnl±^Vn 



1 
2 



2/5 

1-/3 



Sl2 



1-/3, 2(3 



1-/3 



-- In^ + 2 In ^— ^ In ^— ^ + 2 Lia^— ^ + -tt^ 



1 + /3 



1 + /3 1 + /3 



'1 + /3 



2 
3' 



(3.24) 



and 



4iVC, 



1 , 1-/3 /m^ 
— In 



1 + /3 1^12 



(3.25) 



Observe that the sum of the soft and ultra-violet renormalized virtual terms is finite, A\ + 
A\ = 0, as required We are therefore free to return to 4 dimensions with the finite 

remainders Aq + A^. The final result for the 0{as) correction consists of two contributions 
to the cross section: a two-body term a'-^-* and a three-body term cr^^^ where 



(2) 



(3.26) 



and 



0" 



(3) 



1 



J2\M3\'dT,. 



(3.27) 



2Si2 

As a check, these results may be integrated to give a total rate and compared against the 
known next-to- leading order result ||T3] taken from WM: 



Otot 



a. 



ao^Cp (32 - 8p') U2{t) + (16 - 4p^)F3{t) 
47r L 
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FIG. 5: The next-to-leading order contribution to the total cross section for producing a massive 
quark pair in electron-positron annihilation via single photon exchange. The two- and three-body 
contributions together with their sum are shown as a function of the soft cutoff 5s- The bottom 
enlargement shows the sum (open circles) relative to ±5% (dotted lines) of the analytical result 
(solid line) given in Eq. ( 3.28| ). 



+ (2 + p) sjl - p F^it) + (8 - 2p2) ln(t) ln(l + t) 
+ (-12 + 2p + 7pV4) ln(t) + (3 + 9p/2)Jl - p 



(3.28) 

where cxo is the leading order cross section for producing a pair of massless quarks given by 



3Si2 ^ 



-NQt, 



and we have also defined 



with p = Aw? / s\2 and 



F3(t) 

F4(t) 



l + ^T^p 



Li2(-t) + ln(t) ln(l -t) 

6 ln(t) - 8 ln(l - t) - 4 ln(l + t) 



(3.29) 



(3.30) 



(3.31) 



The next-to-leading order corrections are shown in Fig. ^ for -^i =11 GeV and m = rub = 5 
GeV. The two- and three-body contributions together with their sum are shown as a function 
of the soft cutoff 6s- The bottom enlargement shows the sum (open circles) relative to ±5% 
(dotted lines) of the analytical result (solid line) given in Eq. ( ^.28[ ). The result quickly 
converges to the known result. 
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It is satisfying that the fully inclusive rate from the slicing method agrees with that from 
Ref. 



Having made this necessary check, the results may be used to histogram a wide 
variety observables and to study various physics issues. We refrain from any such studies 
here, and instead pass to our second example. 



B. Electron-positron annihilation to massless quark pair 

The process to be studied in this section is similar to that of the last section, but with 
one key difference: the quarks are considered massless from the beginning. Therefore, in 
addition to the final state soft singularities there are final state coUinear singularities. The 
leading order cross section 

da'' = -^Y,\M2\^dT2, (3.32) 

2Si2 

is expressed in terms of the (summed and averaged) matrix element squared 

^|M2p = 2Ne'Ql (^k±^ _ e) , (3.33) 



s 
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calculated from Fig. |T], and the two body phase space 



dT2 = ^( — ] — r / sin^-^^^ de . (3.34) 

16n \Si2/ r(l - e) Jo 

In four space-time dimensions, integration over the phase space produces the result shown 
previously in Eq. ( p.29| ). 

For the QCD corrections we need the soft cross section ( p.5| ), the hard-collinear cross 
section ( |2.25| ), the hard-non-coUinear cross section ( |2.26| ), and the virtual contribution. The 
real emission diagrams are shown in Fig. ^ where the quark lines are to be interpreted as 
massless. The two-to-three body matrix element squared needed to evaluate the hard- 
non-collinear cross section Eq. ( 2.26| ) follows directly from Eq. ( p.5|) of the previous example 
by setting m = and evaluating the traces in four space-time dimensions. The soft cross 
section Eq. ( p^) may also be obtained from the results of the last example by setting m = 
in Eq. (|3.71 ) giving 



dafi = da^ 



as r(l - e) (Ann. 



27rr(l - 2e) V su 



Cf / dS . (3.35) 



The pole needs to be integrated over the soft phase space according to Eq. (|2.17|) . To this 
end, define 

1 



^(■335-545) = / dS . (3.36) 

Using the massless (/3 = 1) form of Eq. ( |3.12| ), the energy integral Eq. ( |3.13| ), and the angular 
integrals given in Appendix B, the result is 

/(S35S45) = f 4 - - In 5, + 2 In^ 6,) . (3.37) 
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We may therefore write the final expression for the soft cross section as 



daQ = da^ 



as r(l-e) /47rp^ 
27rr(l-2e) 1^ S12 



# + - + ^0 



(3.38) 



with 



^2 



2Cf 



(3.39) 



The final state hard collinear cross section was derived in Sec. B 1| . The relevant splitting 
is q —>■ qg. From Eq. ( p. 451) we have 



da' 



rri 



)_ 

27rr(l-2e) 1^ 



At 



+ Ay^' 



(3.40) 



with 



Al 



Cf (3/2 + 2 In 5, 



^0 



7/2 - 7rV3 - In^ 5, - In 5^ (3/2 + 2 In 5, 



(3.41) 



The interference of the one-loop diagrams in Fig. ^ with the leading order diagram yields 
the virtual contribution. In Fig. diagrams (6) and (rf) add to zero via the Ward identity. 
Diagram (c) vanishes for massless quarks. This leaves diagram (a), comprised of the vertex 
shown in Fig. evaluated for massless quarks. The result for the vertex is 



[-ieQq] 5iju{p2)lt,u{pi)-f{q^ 



with 



Zn 



_g2 



r(l -2e) \e 



1 3 , TT 

+ :^ + 4 + 



2e 



6 



(3.42) 



(3.43) 



We may therefore write the final expression for the virtual contribution to the cross section 
as 

'as r(l-e) /47r/i2\ 



day = da 



27rrfl - 2e) 



S12 



A\ 



(3.44) 



with 



Al = -2Cf 
Al = -3Cf 



A"" 
^0 



-2C^(4 - 7rV3) . 



(3.45) 



The full two-body weight is given by the sum das + day + 2da'^'^^ . The factor of two 
occurs since there are two quark legs, either of which can emit a gluon. At this point we 
have a finite result since y4| + = and A\ + A\ + 2Al~^'^^ = as required . The finite 
two-body weight is given by 



a 



(2) 



(3.46) 
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FIG. 6: The next-to-leading order contribution to the total cross section for producing a massless 
quark pair in electron-positron annihilation via single photon exchange. The two-body (negative) 
and three-body (positive) contributions together with their sum are shown as a function of the 
soft cutoff 6s with the collinear cutoff 6c = 5^/300. The bottom enlargement shows the sum (open 
circles) relative to ±5% (dotted lines) of the analytical result (solid line) given in Eq. ( 3.481 ). 
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FIG. 7: The S35 — S45 plane for electron-positron annihilation to massless quarks showing the 
delineation into soft S and collinear C regions. The triangles marked m give vanishing contribution 
for Sc <^ 6s- 
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while the three-body contribution is given by 



a(^) = a^^ = -^/ V|^3prfr3. (3.47) 

A necessary check may be made by integrating these results and comparing with the 
known analytic answer. The contributions from a'-^-* (negative) and a^^^ (positive) and their 
sum are shown in Fig. |^ for ^/s = 500 GeV as a function of the soft cutoff Sg with the 



collinear cutoff 6c = 5^/300. The known result may be found in p6[, for example, and is 
given by 



CT- 



(1) 

tot 



a 



ctoj^SCf, (3.48) 

where o"o is given in Eq. ( ^.29| ). The bottom enlargement shows the sum (open circles) 
relative to ±5% (dotted lines) of the known result (solid line) given in Eq. ( fj.48D . Very good 
agreement is found below 5^ ~ 2 x 10"'^. 

Before proceeding further, it is instructive to examine some issues related to the cutoff 
dependence of this technique. As shown in Fig. |], the answer converges to the known 
result for Sg < 10~^ when 6c = 5<(/300. We have imposed the requirement 6c <^ 6s which 
may be understood by examining the nature of the three-body phase space for this case. 
Neglecting both initial and final state masses, four-momentum conservation yields = 

534 -|- S35 + S45. The soft region is defined by < 5sy^Si2/2 which, taken with Eq. (p.7|), 
can be recast as S45 < 6sSi2 — S35. This is shown as the region S in the plot of S45 versus 

535 in Fig. ^ Two collinear regions defined by the constraints S35 or S45 < 6cSi2 are shown 
as the regions labeled C in Fig. ^. There are two small regions labeled "m" which are 
properly included in the collinear regions C. However, using a fixed upper limit of 1 — 6s 
in calculating the hard collinear contributions (cf. Eq. ( p.35|) ) these regions are excluded. 
They are also not included in the hard-non-coUinear three-body integrations. With some 
effort, it is possible to analytically evaluate the required integrals ( |2.33| ) over the m regions. 
The result (derived in Appendix C) is that occurrences of In (5c In 5^ in Eq. ( p.47| ) are to 
be replaced by In^cln^s — Li2{6c/6s)- From the properties of the dilogarithm function we 
note that the correction term vanishes like 6c/ 6s in the limit of small 6c- Accordingly, the 
contributions from the regions denoted by m in Fig. |^ may be made negligible by requiring 
6c ^ 6s. Of course, as 6c and 6s become smaller the statistical errors on the sum of the 
two- and three-body weights increase. In practice, one must compromise between the errors 
induced by larger cutoffs and the statistical errors. For many calculations it has been found 
that choosing 6c to be 50 - 100 times smaller than 6s is sufficient for answers accurate to 
a few percent. Acceptable ranges for 6s must be determined on a case by case basis, as 
illustrated by the examples shown here. Furthermore, the sign of the deviations as 6s grows 
differs from process to process. 



C. Electron-positron annihilation to photons 

In this section we consider an example fragmentation process, inclusive photon production 
in hadronic final states of electron-positron annihilation, calculated to leading order in the 
electromagnetic coupling a. For pedagogical purposes only the radiation of photons from the 
final state quark or antiquark will be included, i.e., initial state radiation will be neglected. 
This process is different from the previous two examples in that there are final state collinear 
singularities only, and they are removed through the factorization procedure. 
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FIG. 8: Leading order non-perturbative contribution to inclusive photon production via photon 
exchange. 





(a) 



(b) 



FIG. 9: Leading order perturbative contribution to inclusive photon production via photon ex- 
change. 



The diagram for the leading order non-perturbative contribution is shown in Fig. The 



cross section from Eq. ( 2.52 ) is 



dan 



+7X 



E 

q 



dzD^lq{z)doQ 



(3.49) 



To simplify notation it is helpful to write the Born-level total cross section for e^e — > qq 
in terms of that for e^e~ — > /i'^yU" 



We further denote ctq ^ 



(3.50) 



by cr, and ^ ^ by cr^^. Taking into account D^/g{z) = 
D^/g{z), and letting / denote the quark flavor, we arrive at the result for the leading order 
non-perturbative contribution 



1 da 
(Jfif, dz 



2Nj2Q}D,/f{z). 



(3.51) 



Additionally, there are two- and three-body pieces that make up the leading order per- 
turbative contribution. The Feynman diagrams are shown in Fig. Because there are only 
final state coUinear singularities present, the relevant decomposition of the two-to-three con- 
tribution to the cross section is into collinear and non-collinear terms, a = ac + cr-^. The 



collinear term ac is handled as discussed in Sec. |II B 2| . In this case there are no soft singu- 
larities so the soft-collinear terms A^'^ are not present. The two-body piece follows from Eq. 

1 da^"^^ a ~ 

'^N-T.QtD./fi^,Mf), (3.52) 



Cr^^ dz 



/=1 
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once the replacement — a is made. D^/f{z, Mf) is as given in Eq. (in and may be 
expanded using D^/y{x) = 5(1 — x) + (9(a^) and D^/i{x) = 0{a) for i = g,q. The leading 
term in a is therefore 



Jz y 



(3.53) 



where P^'^'^^ is as given in Eq. (|2.65|) with 



Pjg{z,e) 



l + n-zV- ez^ 



(3.54) 



The final result for the two-body piece of the leading order perturbative contribution is 



1 f/fT^^) (y 



a^^ dz 



In 



z{l — z)6cS 



M) 



+ z 



(3.55) 



The complementary non-collinear piece a-^ follows from the matrix element represented 
in Fig. I 



2s Jc 



Ae' 



6 



LT.m\'dT, = —NY.Qjljr 



/ 2g34g £35 S45 
C VS35S45 S45 S35 



(3.56) 



We could now study the cutoff dependence numerically as in the previous examples. How- 
ever, the integration over phase space may be performed analytically and rather straight- 
forwardly so we take this route to demonstrating the Sc independence of the full result. To 
this end, consider, in the virtual photon rest frame, the three-body phase space 



111 d^pi d^ps 
8(27r)5^ E5 



(5(v^ -Es-E^- E,) . 



(3.57) 



Let denote the virtual photon four-momentum. Taking p^ along the z axis and defining 
X4 = IEa^I \fs^ z = 2E^/y/s we may write the four-momenta as 



q = v^(l, 0,0,0) 

P5 



^^(1,0,0,1) 



P4 = X4^-(l, sin6', 0, cos6') . 

Momentum conservation ^=^3+^4+^5 gives 

\/s \/s \/s 

P3 = {E3, sine, 0, -3^4— cos6' - z—) 

The mass-shell condition P3 = can be used to fix E^ as 

E3 = \/xl + z'^ + 2X4Z cos 6 . 



(3.58) 
(3.59) 

(3.60) 
(3.61) 



(3.62) 
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Writing (Pp4 = 2n dcosOEldE^ and d^p^ = AnE^dE^ the phase space delta function may be 
used to perform the cos^ integral. The invariants S34, S35, and S45 may then be written in 
terms of X4 and z 



The three-body piece is now 
= — 

/=1 



o-^M dz 



rif 



S35 = s{l - X4) 
S45 = s{z + X4 - 1) 
S34 = s(l - z) . 



2(1 -z) 



1 — Xa Z + Xa — 1 

\ \ 

— X4){z \ X4 — \) Z + X4—I 1 — X4 



(3.63) 



. (3.64) 



This is to be integrated over the non-collinear C region defined by S45 > ScS and S35 > ScS 
which is equivalent to 1 — z + 6c < x^ < 1 — 6c- The integral may easily be performed. 
Dropping terms of 0{6c), the final result for the three-body piece of the leading order 
perturbative contribution is 



1 da^^^ 



a^^ dz 



27r ^ 



1 + (1 



In 



(3.65) 



Adding the non-perturbative (|3.51|) , perturbative two-body ( |3.55|) , and perturbative 
three-body (|3.65|) contributions we obtain the well known result ||75| , [76| 



1 da 

a^^ dz 



2Nj2Q}\D,/fiz) + ^Q 



zf 



In 



z\l-z)s 



(3.66) 



which is independent of 6c- 



D. Drell-Yan 



Our next example is that of the QCD corrections to lepton pair production in hadron- 
hadron collisions which illustrates the method for handling initial state collinear singularities 



developed in Sec. [II B 3 



The leading order contribution mediated by a virtual photon is shown in Fig. |10|. The 
leading order partonic cross section 



da^ 



2s 



12 



(3.67) 



is expressed in terms of the (summed and averaged) matrix element squared calculated in 
n = 4 — 2e dimensions 



2 ftl^+ ^23 



and the two-body phase space 



dTn 



i2e 



IQtt \S12J r(l — e) Jo 



sm 



l-2e 



ed9- 



(3.68) 
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FIG. 10: Leading order contribution to Drell-Yan production of a lepton pair via photon exchange. 



For the QCD corrections there are five pieces to consider: the finite hard-non-colhnear 
partonic cross section Eq. ( |2.26| ), the mass factorized hard-coUinear cross section Eq. ( p. 741 ), 
which consists of two pieces, the soft part of the initial state factorization counterterms and 
the mass factorization residuals (the G functions), the soft cross section Eq. ( |2.22| ), and the 
virtual corrections. 

Shown in Fig. |TI| are the real emission diagrams that give jMsp. Defining the (summed 
and averaged) matrix element squared as 



'34 



(3.70) 



we find 



^ = tr [izYi^Y] 
X {-tr[^27M(^i 



+2tr [^27" (^2 - - ^^hA ^15^25 

- tr [^27^^ (^2 - - ^5)7.] 7^25} • 



(3.71) 



The hard-non-collinear partonic cross section is obtained by evaluating the traces in four 
space-time dimensions. 

There is a soft singularity when the gluon's energy goes to zero in the real emission 
diagrams. This contributes to the soft cross section presented in Eq. ( p.22|) . The sum runs 
over the initial state quark lines (labeled by 1 and 2). The color connected Born cross section 



12 



-C^rfff". We find 



as r(l - e) / A-Kfi^.^ 



dS . 



(3.72) 



27rr(l-2e) V S12 / 'J \Pi- PhP2- Pb, 
The integration of the pole over the soft phase space measure ( p.l7|) is written in terms of 

1 



/(tl5t 



25 J 



^15^25 



-dS. 



(3.73) 



In the P1P2 center-of-momentum system we take 

\fSV2 



Pi 
P2 



2_ 

s'12 



1,0,. ..,0,1) 



1,0,. ..,0, 



(3.74) 
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FIG. 11: Real emission contribution to Drell-Yan production of a lepton pair via photon exchange. 



Using these together with Eq. ( p.l2| ) we find 



^15 
^25 



-2pi ■ p5 = -^^1^^5(1 - cos 6*1' 

-2p2 ■P5 = -^/s^E5{l + COS 6^1) 



(3.75) 



Using the energy integral Eq. (|3.13| ) and the angular integrals given in Appendix B, we find 

1 



/(ti5t25) = ( 4 - - In 5s + 2 In^ 5, 
2si2 Ve^ e 



(3.76) 



We may therefore write the final expression for the soft cross section as 



dafi = da 



rri 



) /47r/i 

27rr(l -2e) 1^ S12 



2\ 



,2 + , +^0 



(3.77) 



with 



^2 

A' 



2Cf 

-AC pin 6 s 
ACf In^ 5s . 



(3.78) 



Because the quarks are massless, there is a collinear singularity when the gluon becomes 
collinear to either of the initial state quark lines. This singularity is removed through the 
factorization program described in Sec. |11B 3| . The soft- collinear pieces of the initial state 
factorization counterterms are given in Eq. ( [2.74] ) as 



da 



coll 



da^ 



«s rfi 



27rr(l -2e) 



' Aixjil 



Af 



+ ^0 



(3.79) 



with 



A\^ 

ISC 

^0 



(2 In (5, + 3/2) 



Al' 



S12 



(2 In 5, + 3/2) In ^ 



(3.80) 



The one-loop virtual diagrams are shown in Fig. [T^. As in the case for electron-positron 
annihilation to a massless quark pair, diagrams (6) and [d) add to zero via the Ward identity, 
and diagram (c) vanishes for massless quarks. This leaves diagram (a) for which the vertex 
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FIG. 12: Loop and counterterm corrections to Drell-Yan production of a lepton pair via photon 
exchange. 



shown in Fig. ^ is needed, for massless quarks. The result is given in Eq. (|3.42| ). The final 
expression for the virtual cross section is 



da\r = da 



a, F(l - e) [Anfil 



27rF(l - 2e) \ Su 



4^ A'" \ 

,2 + , +^0) 



(3.81) 



with 



Al 



-2Cf 
-3Cf 

-2Cf{4: - 7rV3) 



(3.82) 



At this point we pause to note that the two-body weight is finite: A^ + A2 = and Af + 
Ai + 2Af = 0. The factor of two occurs since there are two quark legs, either of which can 
emit a gluon. 

In addition to the qq initiated processes, there are also qg initiated processes at this order 
of perturbation theory, as shown in Fig. |l^. The singularities are initial state collinear only 
in origin and arise from the Pqg splitting in diagram (6). They are removed by factorization. 
As the Pqg{z) kernel is finite for z = 1 there are no soft singularities. This implies the A^'^ 
terms in Eq. ([2.74]) are not present; only the finite G terms remain. 



The final finite two-body cross section is given by the sum of the residual G terms from 
both the qq and qg initiated processes and the finite two body weights from the qq process, 
-|- -|- 2^0^^. The result, summed over all parton fiavors, is 



a 



(2) ^ ( 

2tt 



J2 / dxAdxsdao \ G f /a{xa, ^'f)Gj /j^^xb, l^f) {Aq + AI + 2Aq) 



+ G f/A{xA, f^f)Gj/j^{xB, fJ'f) + Gf/A{xA,fJ'f)Gj/j^{xB,fJ'f) + {xa ^ Xb) ■ (3.83) 
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(a) 





FIG. 13: Quark-gluon initiated contribution to Drell-Yan production of a lepton pair. 



The G functions are given in Eq. ( p.75|) . The three-body contribution is given by 

^(3)^ ^ J dxAdxBGi/AixA, IJ'f)Gj/BixB, IJ'f)daij , (3.84) 



with the hard-non-colhnear partonic cross section given by 



da 



1 



2si2 Jhc 



_Y,\Mi''^\^dT 



(3.85) 



Physical predictions follow from the sum a^'^^+a^^^ which is cutoff independent for sufficiently 
small cutoffs. The results may be integrated to obtain the total rate for > Qmin 
checked against the known 0{as) corrections 



dQ'^ [ dxA [ dxBl V Gi/A{,XA,lif)Gj/B{,XB,^if)'^j3 (3.86) 
JqVs JqVSxa ii,t^,,- dQ^ 



+ [Gi/A{XA, fJ'f)Gg/B{xB, fJ'f) + Gg/A{XA, l^f)Gi/B{xB, f^f) 



da. 
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i=q,q 



dQ^ 



where is the square of the lepton pair invariant mass and S is the hadron-hadron center of 
mass energy squared which is related to su, the parton-parton center of mass energy squared, 
via Si2 = xaXbS. Defining z = Q^/si2, the 0{as) hard scattering partonic subprocess cross 
sections are given by 



da, 



qq 



a. 



dQ^ 



ao-^GpUil + z' 

ATI 



ln(l - z) 

I- z 



1 - z 



+ 



27r2 



(3.87) 



and 



da, 



qg 



dQ' 



a, 1 f 3 3 r, ^1 , , 

(^o — -\ - + z--z^ + 2P+(z) 
27r2 1 2 2 '^^^ ' 



In ^ - 1 + In 

z fij 



where ctq 



(3.88) 



(3.89) 
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FIG. 14: The next-to-leading order quark-quark contribution to the Dreh-Yan cross section. The 
two-body (negative) and three-body (positive) contributions together with their sum are shown as 
a function of the soft cutoff Ss with the colhnear cutoff dc = Ss/ 50. The bottom enlargement shows 
the sum (open circles) relative to ±5% (dotted lines) of the analytic result (solid line) given in Eq. 



Pgq(z) and Pgg{z) are the splitting kernels given in Eqs. (|2.55| ) and ( |2.58| ). 

We show numerical results on the cutoff (in) dependence for proton-proton collisions at 
= 28.28 GeV with Qmin = 10 GeV. Hard scales are set to the lepton pair mass /i/ = 
Hr = Q and the number of flavors taken to be = 3. 

Shown in Fig. H is the next-to-leading order quark-quark contribution to the Drell-Yan 
cross section. The two- and three-body contributions to the cross section (negative and 
positive, respectively), and their sum are shown as a function of the soft cutoff 6s . The 
coUinear cutoff 6c = 6s/50. The bottom portion of the figure shows the sum (open circles) 
relative to ±5% (dotted lines) of the analytic result (solid line) given in Eq. (|3.87|) . 

Finally, we show the next-to-leading order quark-gluon contribution to the Drell-Yan 
cross section in Fig. The two- and three-body contributions and their sum are shown as 
a function of the coUinear cutoff 6c- The bottom enlargement shows the sum (open circles) 
relative to ±5% (dotted lines) of the analytic result (solid line) given in Eq. (|3.88|) . 

In both cases, nice agreement is seen with the known analytic result, providing a cross 
check on the use of the two cutoff phase space slicing method. 



E. Single Particle Inclusive Cross Section 

Our final example is that of the single particle inclusive cross section in hadron-hadron 
collisions. The input needed for this calculation includes the squared matrix elements for 
the 2^3 subprocesses f?^ and the results for the 0{al) one-loop contributions to the 
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FIG. 15: The next-to-leading order quark-gluon contribution to the Drell-Yan cross section. The 
two-body (negative) and three-body (positive) contributions together with their sum are shown as 
a function of the colhnear cutoff 5c- The bottom enlargement shows the sum (open circles) relative 
to ±5% (dotted lines) of the analytic result (solid line) given in Eq. ( |3.88| ). 



2^2 subprocesses |78|. For the purpose of this example, the notation of |^ will be 
used, since much of the input needed can be found in the appendices of that paper. The 
partons are labelled + + 2 and A + B ^ I + 2 + 2, for the 2 ^ 2 and 2 3 

subprocesses, respectively. A flavor label is used to denote the flavor of parton A, and 
similarly for the other partons. 

The lowest-order contribution to the inclusive cross section for producing a hadron h in 
a collision of hadrons of types A and B can be written as 



da 



B 



1 



2xaXbs 



w{aA)w{aB) 



E 



GaA/A{xA)GaB/B{xB)Dh/ai{zi)dXA dxsdzi 



(4) 



d,'p)dT2 



(3.90) 



where a = {0^,0^,01,02} and p = {Pa^Pb,PuP2} denote the sets of flavor indices and 
parton four- vectors, respectively. The factors appearing in the spin/color averaging are 
given by 

. . f 2(1 - e)V a=gluon 

I 2N a=quark or antiquark 



wia] 



with = 3 and V = N"^ — 1. The factor dr2 is the differential two-body phase space element 
from Eq. (|2.11| ). Eq. (|3.90| ) gives the contribution where parton 1 fragments into the hadron 
h. Care must be taken to explicitly include in the sum over a those terms corresponding to 
the case where parton 2 fragments into h. For compactness, these terms will not be explicitly 
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written. The squared matrix elements for the various subprocesses, denoted by ip^^^djp), 
may be found in Ref. . 

Next, consider the one-loop virtual corrections to the 2 — 2 subprocesses. These take 
the form 



da" 



7; GaA/A{xA)GaB/B{xB)Dh/aAzi)dXAdXBdZi 

ZXaXbS aA,aB,ai,a2 



X 



{Anas 



w{aA)w{aB) 



a. 



47r/i' 



rfi-e) 



2iT \2pa-PbJ r(l-2e) 



tp'^'\d,p)dr2 



(3.91) 



where 



1 



,PA ■ PB 



+ - y In 
2f ^ 

m,rt 

- ^T.^^'\d,p)+^ljPsia,p} + Oie). 



(3.92) 



This expression for ip^^^ differs slightly from Eq. (35) in Ref. |^ because we have chosen 
to extract a different e dependent overall factor: a factor of r(l + e)r(l — e) ~ 1 + 
has been absorbed into the above expression for ijj^^K Furthermore, the arbitrary scale Q'es 
used in Ref. has, in order to match the conventions used elsewhere in the present work, 

been chosen to be 2pA -pB- The expressions for the functions ip^n V'/^s be found 
in Appendix B of Ref. [0. The quantities C(a„) and 7(0^) are given by 



C{a) 



and 



7(a) 



= 3 a=gluon 
Cp = ^ a=quark or antiquark 



(lliV-2n/)/6 a=gluon 

3Cf/2 a=quark or antiquark 



It will be convenient for subsequent expressions to adopt the following notation: 



2pa-PbJ r(l-2e)' 
The one loop virtual contributions can now be written as 

1 



(3.93) 



da" 



^aA/A (xa) Gag /B{xB)Dh/aAzi)d^A dXB dzi 



X 



2xaXbs 



a A, a. 3,11, 12 
\2 



wiaAjvoiaB] 2t[ \ e ' 



(3.94) 



where 



A\ = - Y.C{anW\a,p) 



(3.95) 
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Al = -E7(a„)^(^)(a,P) + ^ E 1^ (P^) ^p^^fJid^p) (3.96) 

„ Z m,n \PA ■ PB ) 

(3.97) 



^0 = --Y.CM^^'\a,p) +ijfsia,p). 

^ n 

Next, the contributions from the 2 — > 3 subprocesses in the hmit where one of the final 
state gluons becomes soft are needed. The contributions of the 2 — 3 subprocesses may be 
written as 



da 



2^3 



{Anas 



2xaXbs w(aA)w(aB 



-T E GaA/A{xA)Ga/B/B{xB)Dh/ai{zi)dxAdXBdZi 



a A, a. 3,0.1,0.2, 0-3 

X '^{aA,aB,aua2,as,p^,p'^,p1,pi^,p'^)dT3. 



(3.98) 



The expressions for the 2 — >■ 3 squared matrix elements appearing in Eq. (|3.98| ) may be 
found in Ref. ||7S|. As noted earlier for the two-body contributions, one must include in the 
sum all possible parton to hadron fragmentations. 

Consider the case where the soft gluon is parton 3. In this limit, the function may be 
expanded as: 



^{aA,aB,aua2,a3,p%p%,p'(,pf^,p^) ~ E ^os,^ 



Pm ' Pn 



Pm ■ PsPn ■ P3 



-V'm^ (oA,aB,ai,a2,PA,PB,Pi,P2) 



(3.99) 

Next, one must integrate over the soft region of phase space defined by < 
5s\/2pA ■ PbI'2- This is easily done using the integrals given in the appendix. The resulting 
soft contribution may be written as 
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X 



2xa'3^bs 
(47ras 



E Ia{xa) Gas /B{xB)Dh/ai{zi) dXA dXB dzi 



aA,aB,ai,a2 
\2 



as Ml , Al 



w{aA)w{aB) 27r V e 



+ ^ + A'o dV 



where 
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J2G{a^)^^'\d,p) 

n 

-2in5,E^K)^^'H«,P) - ^ E In f^^l i^l^fJ{s,P} 

n ^ \PA-PbJ 

21n2 5,EC(an)^^'Ha,P) 

' P2 ■ P3 ' 



(3.100) 

(3.101) 
(3.102) 



(4,c) , {4,c) 



+ «2^+<^) 



iinM^i:^ 

2 \Pa-Pb, 
lln^ (P^^ 



.Pa-Pb, 



+ Li2 
+ Lia 



^Pa ■ Pb , 
' P\-p-i^ 
J>a-Pb, 



+ 21n5^ In 
+ 21n5o In 



P\ -Ps 

,Pa-Pb, 

' P2-P3^ 

,Pa-Pb, 



(.3.103) 



After the collinear singularities associated with the two parton distribution functions 
and the fragmentation function have been factorized and absorbed into the corresponding 
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bare functions, there will be soft-collinear terms left over due to the mismatch between the 
integration limits of the collinear singularity terms and the factorization counterterms. In 
addition, there can be collinear singularities associated with the non-fragmenting parton in 
the final state, corresponding to gluon emission or qq production. Collecting together both 
types of collinear terms, the result can be written as follows: 

da""^^ = — ^ GaA/A{xA)GaB/B{xB)Dh/ai{zi)dXAdXBdZi 

ZXaXbS aA,aB,ai,a2 

.^i-)l^^|i(^,<-),r. ,3.104) 

where 

Af' = 5][21n5,C(a„)+7(an)] (3.105) 

n 

A™" = ^[21n5,C(a„)+7K)]lnf^^^) +[21n5,C(ai)+7(ai)]lnf^^f-^ 

A,B \ J \ / 

+ i\a2). (3.106) 

Here /x/ and Mf are the initial and final state factorization scales. The function 7'(a) is 
given in terms of the hard collinear factors of Eqs. ( |2.46| ) - ( |2.51| ) as 

s I A'K^''^ a= quark or antiquark 

^ = I a™ + Ar"^" a = gluon ^3.107) 

After the mass factorization has been performed, the bare parton distribution functions 
and fragmentation functions have been replaced by scale dependent MS functions. In addi- 
tion, there are finite remainders involving the G and D functions: 

da = - — ^ V (47ras) — ^(^x^dxBdziip'-^\a,p)dr2 

2xaXbS aAm,a2 w{aA)w{aB) 271 



X 



Ga^/AixA, IJ-f)GaB/BixB, fif)Dh/a,iZl, M f) 
+ GaA/A{xA, ^^^f)Gag/B{xB, ^J^))Dh/aAZl^ ^/) 
+ GaA/A{xA,^J^f)GaB/B{xB,^^))Dh/aAZl^^})\■ (3.108) 

At this point, all of the singular terms have been isolated as poles in e or have been 
factorized and absorbed into the bare parton distribution and fragmentation functions. The 
e dependent pole terms all cancel amongst each other: 

Al + Al = Q (3.109) 
Al + Al + = 0. (3.110) 

The finite two-body contribution is given by 

da"^^"^ = da^ + da 



X — 
27r 



Al + Al + Af^ dxAdxBdzidT2. (3.111) 
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The three-body contribution, now evaluated in four dimensions, was given in Eq. (p.98|) 
where now the soft and colhnear regions of phase space are excluded. 

As in the previous examples, the structure of the final result is two finite contributions, 
both of which depend on the soft and colhnear cutoffs - one explicitly and one through the 
boundaries imposed on the three-body phase space. However, when both contributions are 
added while calculating an observable quantity, all dependence on the cutoffs cancels when 
sufficiantly small values of the cutoffs are used. 

IV. DISCUSSION AND CONCLUSIONS 

A technique for performing next-to-leading-logarithm calculations using Monte Carlo 
techniques was described in detail. The method uses two cutoff parameters which serve 
to separate the regions of phase space containing the soft and colhnear singularities from 
the non-singular regions. The main derivations for experimentally degenerate, tagged, and 
heavy quark final states were given, as was a discussion of initial state factorization. We 
provided five illustrative examples applying the method. 

The first example was that of the QCD corrections to electron-positron annihilation into 
a massive quark pair. The quark mass serves to regulate any would-be final state colhnear 
singularities. The final state soft singular region is delineated using one cutoff. The second 
example was that of QCD corrections to electron-positron annihilation into a massless quark 
pair. In this case final state soft and colhnear singularities are encountered. Both soft and 
colhnear cutoffs are therefore required. They should be chosen such that 5c -C 5s- The 
third example of inclusive photon production in hadronic final states of electron-positron 
annihilation was presented, illustrating the use of fragmentation functions. Finally, the QCD 
corrections to lepton pair and single particle production in hadron-hadron collisions were 
given. These examples include both initial and final state soft and colhnear singularities. 
The use of scale dependent parton distribution and fragmentation functions was explained. 

The Monte Carlo results, integrated to give an inclusive cross section, were shown to 
be in complete agreement with those available in the literature. This is not the end of 
the utility of the method, but only the beginning. Given the full access to the parton four 
vectors and corresponding weights, we are free to combine them in any way that is consistent 
with an infrared-safe measurement function, which may include a jet finding algorithm and 
experimental cuts. 

The method has been applied to a wide range of hard scattering processes and it has 
been found to be both simple to implement and numerically robust. 
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APPENDIX A: COMPARISON WITH OTHER METHODS 



The essential difference between the phase space shcing and the subtraction methods may 
be gleaned from the following simple example l^^. Consider an integral to be calculated: 




- -F(0)| , (Al) 
ox e J 

where F{x) is a known but complicated function related to a two-to-three body matrix 
element. The variable x represents either the energy of an emitted gluon or the angle 
between two massless partons. In a traditional fully or single particle inclusive calculation 
the integral / would be performed completely analytically. 

In the subtraction method one simply adds and subtracts F{0) under the integral sign. 

/ = lim I —x' [F{x) - F(0) + F(0)] - -F(0)| 

e^0+ [Jo X e J 

= /' - - F(0)] , (A2) 

JO X 

giving a finite and numerically calculable result. No approximations are made, however 
in any numerical implementation there will necessarily be a lower limit related to machine 
precision below which the integral must be cutoff. This is not a problem in practice. 

In the phase space slicing method, the integration region is divided into two parts < 
X < 6 and 6 < x < 1 with 5 -C 1. A Maclaurin expansion of F{x) yields 

—x'F(x) + C —x'F(x) - 
ox Js X 

+ F(0)ln5 + C(5). (A3) 

Clearly, the parameter 5 must be chosen small enough so that the term linear in 5 may be 
neglected. At the same time it must not be so small as to spoil the numerical convergence 
of the first term. 





APPENDIX B: SOFT INTEGRALS 



In evaluating the soft integrals we encounter angular integrals which may be written in 
the form 

/(^'') = r d9, sin"-3 9, r dO, sin"-^ 9, ^^ J'^V'Z^'K ^ • (^l) 

" Jo Jo {A + Bcos9i + Csm9iCos92y 



A large collection of these appear in the appendix of Others may be found in the 

appendix of ||4^ or else computed as explained in ||70|]. Here we collect together the results 
covering most of the cases encountered using the two cutoff slicing method. The first two 



are from |4|] with + 



j{o,i) 



TT 



VB^ + C2 
-(n-4) 



In 



Li2 



'A+VB^ + C^ 



A-VB^ + C^ 



fB2) 



+ - In 



' A + ^B^ + C^^ 



A+VW+C^J 4 \A-VW+C^ 
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/(0,2) 



27r 



1 - 4) , = In \ 



where we drop 0((n— 4)^) terms. The second two are from [69| with b = —a. If A'^ 

1 1 /A + 5\"/2-2r 



aAn-4 
whereas if ^ B"^ + 



2A 



TT 



a(A + S) [n-4 



+ In 



. 1. .n2t. fA-B 



{A + Bf 



'A-VB^ + C^^ 



A + B 



'A + VB^ + C^' 



A-VWTc'^ 



2 Lis 



B + vWTc^^ 



A-VW+C^ 



2 Li. 



'B-VB^ + C^^ 
A + B 



(B3) 
B^+C^ 

(B4) 



(B5) 



again dropping 0{{n — 4)^) terms in the second of these. The dilogarithm function Li2(a;) 
is defined in []7D[ and numerous useful properties are summarized in . 



APPENDIX C: RECOVERING THE 0{6c/6s) TERMS 

In this appendix we integrate the Pqq{z, e) splitting kernel over the hard-coUinear portion 
of phase space for the case of a 45 singularity as it pertains to the discussion given at the 
end of Sec. |1II B|. Recall that this region is defined by 



hard : Sg— — — < < 

2 - ^ - 2 

coUinear : < S45 < ScSu . 



(CI) 



From Eq. (|2.29| ) we have S34 = (ps + p^y = 2p3 ■ ~ (2p3 ■ ^45)2; and = (ps + ^45)^ - 

S45 + 2p3 -^45 which together yield S34 ^ z{si2 — S45). Using = (si2 — S34)/2y'si2 the hard 
condition becomes 

0<z< . . (C2) 



S45/ s 



12 



The approximation made in Sec. P^I B 1| was to set S45 = in the denominator, in light of the 
collinear condition. This resulted in a decoupling of the z and S45 integration limits in Eq. 
( p.36| ). Relaxing the S45 = approximation gives rise to terms 0{6c/Ss) as now described. 
Keeping the S45 dependence, the required integral is 



•^-^la dS45 /S45 



dz[z{l- z)]-'P,,{z,e). 



S45 \Sl2' 

We may expand Pgq{z, e) about e = and make a change of variables u = S45/S12 giving 

rSc 



(C3) 



I = Cp I u 
10 



-l-e 



F(u) du . 



(C4) 
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with 



F(u) 



F{u) may be evaluated with the help of 



dz 



1 + 



a 1 + 



pa 

/ dz{l-z) 
Jo 



* \nz 

dz 

1 — z 

ln(l - z) 
1-z 

z^hiz 

dz 

1 - z 

z'^ ln(l - z) 



a 1 - 



dz 



Li2(l - a) - — 
o 

-iln'ri - a) 
2 ^ ^ 



e(l -z) -eln [^(1 - 



21n(l - a) 
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1 + z' 



TT 



a lH — I — a lH — )lna + Li2(l — a) 



dz 



2 

ln(l — a 



- - a 1 + - 



6 

I ln(l - a) 



The resulting terms in F{u) may be integrated over u using 

- + In 5^ 



duu ^ " ( 

1-6. 



0,1,2 



-l-e 



In 1 







1-u 



h In 4 ) In (5s - U2{5c/5s) 



for the terms multiplied by 0{e^) in Eq. (|C5|) and 

rSc 

Jo 



duu-^-'lii^ 1 



-1-e 



1-^; 

1 -M 



duu ' 



1 - n 



In I 1 - 
In 



duu 



-l-e 



Li2 1 



1 




1 


— u 


1 




1 


— u 


1 




1 


— u 


1 




1 


— u 



e 



-- In 6s; i = 1,2 
e 



0; 




1,2 



(C5) 

(C6) 
(C7) 
(C8) 
(C9) 
(CIO) 
(Cll) 

(C12) 

(C13) 



(C14) 
(C15) 
(C16) 
(C17) 



for the terms multiplied by 0{e^) in Eq. (|C5|) . Terms containing or leading to contributions 
of 0{6c) or 0{6s) have been dropped. Taking the coefficients of 1/e and e° gives the desired 
result 



(3/2 + 2 In 5, 



7/2 - 7rV3 - In^ 5, + 2 U2{6c/Ss) - In 4 (3/2 + 2 In 5,) 



(C18) 
(C19) 



The second equation is identically Eq. ( p. 471) with the addition of the advertised Li2{6c/6s) 
term. A similar analysis may be performed for the Pgg splitting case with the same result: 
In 6c In 6s In 6c In 6s — Li2 {6c/ 6s ) . 
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APPENDIX D: IMPROVING CONVERGENCE OF TILDE TERMS 



We want to demonstrate how the numerical convergence of the D and G functions may 
be improved. To this end consider the integral 



I —Gj{x/y,fi)Pij{y) , 

Jx V 



with 



ix y 
Pij{y)^Pij{y)\n(s, 



(Dl) 



(D2) 



y /^/ 

Here a logarithm of Sg is numerically being built up. Convergence will be improved if we 
rewrite the result in a form where the logarithmic dependence on 5s is manifest. To do so, 
use the fact that 

l]m[{l-y)PM]^2CAj, (D3) 

y-*l 

with Cg — N and Cq — Cp- Now add and subtract the leading singular piece under the 
integral sign: 



dy 

X y 



Gj{x/y,fi)Pij{y) In 5, 



1-ys 



12 



y /^/. 



-G,ix/y,^i)P:Ay) 



i-ysi2 



Gj{x,ij)-^\n Sc + Gj{x,ij)-^\n 5. 

1-y \ y i^f I 1-y V 



1 - y gi2 
y /^/. 



(D4) 



Regrouping terms gives 
'■i-'^s dy 
Ix y 



Jx 



Gj{x/y,iJ,)Pij{y) - Gj{x,fi)- 



2Ci6ij 



In U, 



2Ci6i. 



Gj{x/y,t,)Pr{y)+Gj{x,^)j^\n \^6, 



The last term may be evaluated with the help of 



l-5s 



dy / 1-1/' 

in a 



- y) 

The final desired expression is 
Fi = Ci5ijGj{x,n) 
dy 



1,2/ 
— m ' 

2 



y 

1 - y si2 \ 

y /^// 

l-Sj 2 



1 - 1/ -§12 

y /^/. 



a(l — x) 



X 



(D5) 



(D6) 



In^ S, 



c 2 
llj X 



+ / —{ Gj{x/y,n)Pij{y)-Gj{x,n) — 
Jx y \ i 



In 5. 



1 ~ y gl2 

y /^/. 



(D7) 



The In 5s is now evident in the first term, and absent from the second integral term. Nu- 
merical convergence will therefore be greatly improved. 
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